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I. 


Introduction and Research Objectives 


An Important class of turbomachinery noise Is produced by the 
Interaction of the blades of a machine with flow perturbations within 
the machine. Such flow perturbations may be stationary in time aitd 
position, as would result from the potential field of a downstream 
strut. Disturbances originating from struts exist continually during 
operatiou of a machine, and will contribute to tone noise depending on 
the strength of the consequent flow response and pressure fluctuations. 

The VPI&SU Turbomachinery Research Group had for some years carried 
on an investigation which was designed to provide insight into the 
fundamental aspects of fan rotor-downstream strut interaction (1,2,3). 
High response, miniature pressure transducers were embedded in the rotor 
blades of an experimental fan rig. Five downstream struts were placed 
at several downstream locations in the discharge flow annulus of the 
single-stage machine (Fig. 1), Significant interaction of the rotor 
blade surface pressures with the flow disturbance produced by the down- 
stream struts was measured. Several numerical procedures for calculat- 
ing the quasi-steady rotor response due to downstream flow obstructions 
were developed (3). A preliminary comparison of experimental and calcu- 
lated fluctuating blade pressures on the rotor blades shows general 
agreement between the experimental and calculated values. Although 
progress has been made as a result of this effort, there are still areas 
in which further work is needed. 
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The objective of the program Is to investigate the interaction of 
fan rotor flow with downstream struts. The results of this investiga- 
tion will be an improved understanding of noise prediction with the goal 
of developing a design-for-noise capability and development of a data 
base for use in validating analytical prediction methods and noise- 
reduction concepts. The experimental program will focus on further data 
reduction of the existing data base from previous work. 

The signal from the blade mounted pressure transducers will be used 
to calculate the unsteady lift and moment coefficients. In addition, 
power spectral analysis on the time-averaged data will be performed to 
retain the phase Information. Both of these analyses will be used to 
validate the tisae-marching code, which was written to obtain the rotor 
blade flow response to a downstream strut row. 

The present report gives a summary of the work for the period 
6/16/84 - 3/15/85. Major topics to be covered in this report include 
the data transfer from NASA Langley Research Center to VPI&SU, further 
data reduction and the calculation of unsteady lift and moment. The 
report will conclude with a discussion of future plans. 
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II. Data Transfer NASA Langley Research Center to VPI&Sti 


This section describes in detail the transfer of analog data re- 
corded using FM tape recorder to digital data that can be reduced using 
computers at VPI&SU. A listing of the programs Is provided in the 
Appendix. 

I. Digitization of Analog Data 

The original analog data were recorded on half-inch tape running at 
a speed of 15 inches per second (ips). The data were digitized in the 
Central Data Transcription Facility at NASA Langley Research Center, 
Analysis and Computation Division, Computer Management Branch, Data 
Management Section. The data were digitized with an anti-aliasing 
filter set at 5 kHz flat. The digitization gave 360 points per cycle of 
rotor revolution. Since the rotor rotational speed fluctuates with 

time, the corresponding data sample rate will have to be adjusted to 
account for thi 3 . The time increment between data point varies between 
50 and 70 microseconds. The once per revolution signal is used to 
determine the sampling rate to give 360 points per rotation. The digi- 
tized data were written on magnetic tape using the UNIVAX VARIAN Model 
77-600 computer. The specifications on the format of the data written 
on the magnetic tape include the following: Nine tracks; Non- labeled; 

Edcdic; 1600 BPI; Fixed block binary; 16 Bit Integer Words and Two's 
compliment (Integer *2). The characteristics of the digitized data file 
on magnetic tape include the following: Blocksize 3840, Record length 

3840 and Record format F. 
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2. Transfer of digitized data to computer at VPI&SU 

The digitized data is transferred to an IBM 3081 system. The 
system operates under the control of VM (Virtual Machine) and OS/VS 2.3 
(MVS) and has twenty-four (24) megabytes of memory. A series of soft- 
ware programs were developed to facilitate the data transfer. Figure 2 
shows the flow chart to read the data off the magnetic tape. The first 
program, TPP CNTL, executes a system program (TPPRINT) which checks the 
tape characteristics. The program TPPRINT scans a tape and tells the 
number of files, the block storage size, and the number of records per 
block (thus record length). Record length is equal to blocks lze/ (number 
of records per block). TPPRINT also gives a machine language transla- 
tion of the stored data, which is not readily readable. 

Next step in the data transfer is to copy NASA's tape to a VPI 
system tape. The VPI system tape is permanently active, easier to 
access and is in readable format. Program TAPE TAPE is written to 
create the system tape. 

Two programs were written to read from the system tape. TAPE L1STE 
shows the first 20 elements of every line. This is not a necessary 
operation, but allows a visual look at a swath of the whole file from 
top to bottom. The second program, TAPE LISTL, shows the first few 
lines of the file in totality. This allows viewing the arrangement of 
data on each line. Some prior knowledge of the data expected aids in 
selecting the format used in TAPE LISTL. (See appendix for a sample 
output from the program TAPE LISTL). 

As mentioned earlier, the time Increment between data point varies 
between 50 and 70 microseconds. The rotor test speed of 2900 RPM flue- 


tuates. The analog data were digitized by 1) assigning a square wave 
pulse to match the once-per-rev nignal, 2) dividing each pulse into 360 
units i and 3) making a digital reading of each channel at each of 360 
cycle positions in every revolution. Thus ATiroe = APosl cion/ rotor 
speed. Since the rotor speed fluctuates, the time increment fluctuates 
and thus cannot correlate exactly with position. It is decided that the 
data will be analyzed in the position domain. 

Next, the square wave must be checked to see that 360 readings 
occur In each cycle. The file POS CHECK searches for the drop In each 
square wave (6000 + followed by 600*) and counts the number of entries 
since the praceedlng square wave. It then shows the number of cycle 
positions In each revolutions, some truncation may occur and some cycles 
show 359 or 361 positions Instead of 360. The approach taken here is to 
discard the bad cycles and the program CYCLE SELECT Is written to do 
that. The program also Ignores the first Incomplete cycle and provides 
a printout of data points retained. A similar file, VOLT PLOT, performs 
the same data manipulation and then plots each channel, the square wave, 
and the once-per-rev signal versus cycle position for three good cycles. 



















III. Other Accomplishments 

A. Presentation and Publication 

"Calculation of Unsteady Fan Rotor Response Caused by Down- 
stream Flow Distortions,” by O'Brien, W. F., Jr., Richardson, S. 
M., and Ng, W. F., AIAA Paper No. 84-2282, presented by W. F. Ng at 
the AIAA 9th Aaroacoustlcs Conference, Williamsburg, Virginia, 
October 15-17, 1984, submitted for publication to AIAA Journal of 
Propulsion and Power . 

B. Improvement In Analytical Code 

Refinement In the time-marching code has led to Improvement In 
the convergence of the solution. The calculation was carried out 
to 800 cycles with error In mass conservation less than 0.7% and 
error In total pressure conservation within 0.2%. The predicted 
pressure fluctuations on the blade surfaces compared very well with 
the measured data [4]. In addition, the Douglas-Newman program and 
the time-marching code are now better documented and readily avail- 
able for future use. 

C. Data Reduction 

All analog data recorded on tapes had been digitized and 
transferred to the computers at VPI&SU for further data reduc- 
tion. A computer program was written to take ensemble average on 
the data. Figure 3 shows some typical raw data and 200 cycles 
ensemble averaged data from pressure transducers mounted on the 
pressure and suction sides of the blade, with and without the 
downstream struts. The rotor Is running at 2900 RPM and 15% surge 
margin. The transducers are at 50% span. 
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A progress 1b also developed co calculate the unsteady lift and 
moment from the blade mounted transducers. The lift Is calculated 
using 


l ■ X <p p ■ p »’‘ m ‘ <i> 

where P n and P. are the pressures on the pressure and suction 
p s 

surfaces of the blade, AA Is the differential area covered by each 
transducer, the pressures can be written as two components: 


and 


V '> ■ P P + p p U) 

p s <t) -?; + p^t> 


where T is the steady-state pressure and P' is the time-varying 
pressure. Thus 


1*1 


<v 


Vi 


M t + 


i»i 


(P . . P .) 

p si 


AA, 


( 2 ) 


The first term is a constant and does not change with time. 
The second term can be calculated from the measurement using the 
blade mounted transducers. Figure 4. A shows the unsteady lift 
against rotor position for different downstream strut locations. 
Similarly the unsteady moment is given by 
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M - l d.(P n - P )AA + I d. (P* - Pi)AA (3) 

1-1 1 p 8 1 1-1 1 P 3 1 

where is the distance between the transducer and quarter 
chord from the leading edge of the blade. 

Again the first term does not change with time and the second 
term can be calculated from the blade mounted transducers. Figure 
4-B shows the unsteady moment plotted against rotor position for 
different downstream strut locations. 

D. Interactions with Industry 

There is strong interest on the current work from the Aircraft 
Engine Business Group of General Electric Company at Cincinnati, 
Ohio. Dtscuasions with their aeroacousttc design group have Indi- 
cated that the problem is still of definite interest to them and to 
the aircraft manufacturer. They may use the analytical code for 
future investigations. Two specific areas of Interest arc: (1) 

Optimum positioning of the stators between the rotors and the 
struts to minimize unsteady Interaction, (2) Effect of nonuniform 
circumferential spacing of the support struts on the fan-strut 
Interaction (5). There will be continuous interactions between the 
Turbomachinery Research Group and aircraft engine Industry. 












Pressure (Pa) 



Fig. 3.B.2. Unsteady pressure on suction surface of rotor (raw data) with strut removed 
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Pressure (Pa) 



Unsteady pressure on pressure surface of rotor (200 ensemble averaged) with strut removed 










IV, Summary and Future Plan 

All analog data recorded on tapes have been digitized and success- 
fully transferred to the computers at VPI&SU for further data reduc- 
tion. In addition, a program was developed to calculate the unsteady 
lift and moment from the blade mounted transducers. The result of this 
analysis will be used to validate the numerical calculation. Work is in 
progress for implementing a digital Fast Fourier Transform (FFT) and 
Discrete Fourier Transform (DFT) subroutines for looking at spectral 

analysis of the data. A parametric study on the effect of different 
Fourier Transforms on the spectral analysis is underway. 

In order to isolate the effect of downstream struts from other 
harmonic disturbances, it is necessary to subtract the amplitude of 
pressure fluctuation due to "background" disturbances from the amplitude 
of pressure fluctuation due to’ the .downstream struts alone. In the 
previous data analysis, this was done by performing a power spectral 

analysis (FFT) on the raw data, then taking 200 averages of the power 
spectrum to give an overall spectral analysis. This was done for both 
cases; with and without the downstream struts. The difference in ampli- 
tude of the pressure fluctuation was then reported in Ref. (1). No 

phase information was retained in this data analysis. 

Further data reduction will include an enkanced-signal spectral 
analysis (6) for the data set. The raw data will first be ensemble 

averaged (conditional averaging) and then a power spectral analysis will 
be performed on the averaged data while retaining the phase informa- 
tion. The vector difference in amplitude of the presure fluctuation 
(with and without the downstream struts) can be obtained using a polar 


plot (Nyquist plot). The result of this analysis will be compared to 
previous results as well as to the analytical model. 

In addition, using the existing analytical code, an effort will be 
Initiated to Investigate the optimum positioning of the stators between 
the rotors and the struts to minimize unsteady Interaction. 
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[ '] Appendix: Listing of programs used In data transfer 

1. Program XPP CNTL - execute a system program (TPPRINT) which checks 
i ' the tape characteristics. 


j 



TPP CNTL 


//B0589TLO JOB 529E5.*TIM OLSEN* , REGION® 1Q24K, TIME=< # 59 ) 
/■JOBPARM LINES»99 
/■PRIORITY STANDARD 
/■LONGKEY WINOFAL 
//STEP0001 EXEC PGM a TPPRINT 
DD SYSOUT®* 

DD 


//SYSPRINT 

//SYSUT1 

// 

// 

// 

// 

// 

//SYSUTZ 
//SYSIN 


DD 

DD 


DSN®VPI13, 
DISP s 0LD« 
UNIT a TAPE# 
LABEL a (I«NL)» 
V0l a SER=VPI13, 
DCB a BLKSIZE® 32000 
SYSQUT®* 


SYSTEM PROGRAM WHICH CHECKS TAPE 
CHARACTERISTICS. 


FILES a 99 1 STOPAFT «3 


LARGEST POSSIBLE; TAPE KILL PROBABLY 
HAVE SMALLER BLOCKSIZE. DSN AND VOL 
GIVE NAME OF NASA TAPE BEING CHECKED. 


/* 

// 


u 

0 

0 

0 
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. Program TAPE TAPE - transfer NASA's tape to VPI’s system tapa. 


TAPE TAPE 


'DATA* DIMENSION IS MALE OF 
THE TAPE'S RECORD LENGTH. 
CALL PL 1 COMPILER. 

PL 1 PROGRAM. 


) ; 


//B0589TLO JOB 529E3.'TIM OLSEN' , REGION" 1536K> TIME a < > 59 ) 

/■LONGKFY HINGFAL 
/■ROUTE PRINT VPIVM1 . OLSENTI. 

/•ROUTE PUNCH VPI VM1 . OLScNTL 
/ ■ JOBPARM L I N E S * I 0 
/■PRIORITY IDLE 
//STEP0001 cXEC PLIXCLO 
//PLI.SYSIN DD ■ 

■ PROCESS SOURCE. FLAO(I). GOSTMT ; 

READTP . 

PROCEDURE OPTIONS! MAIN ) ; 

DECLARE SYSPRINT OUTPUT PRINT FILE ; 

DECLARE 5YSIN INPUT RECORD FILE ; 

DECLARE OUTFIL OUTPUT RECORD FILE ; 

DECLARE DATA! 1920 ) FIXED BINARY! 13 ) ; 

DECLARE COUNT FIXED BINARY! 31 ) ; 

DECLARE END_OF FILE FIXED BINARY! 31 ) INITIAL! 0 
ON ENUFILE! - SY5lN ) END OF FILE * 1 ; 

ON ERROR SNAP BEGIN ; 

ON ERROR SYSTEM ; 

PUT DATA ; 

STOP ; 

END ; 

DO WHILE! END OF FILE * 0 ) ; 

READ FILE! SY5lN ) INTO! DATA ) ; 

WRITE FILE! OUTFIL ) FROM! DATA ) ; 

END ; 

END READTP j 

/« VOL SHOWS OLD TAPE NAME; LABEL 8 DSN SHOW FILE'S POSITION ON TAPE. 
DD D5N=F1LE1 , 

DISP*OL D, 

UN IT “TAPE. 

LABEL*! 1 .NL ) . 

VOL*SER*VPI17 , 

DCB = !RECFM = F»LRECL*38A0.BLKSIZE' , 3840) 

//■ SYSIN REFERS TO OLD NASA TAPE; OUTFIL REFERS TO NEW SYSTEM TAPE 
//* DSN SHOWS ASSIGNED FILE NAME ON NEW TAPE; LABEL SHOWS FILE'S 
//■ POSITION ON TAPE; SER= ' SERIAL • ASSIGNED BY SYSTEM FIRST TIME 
//« AND SHOWN IN OUTPUT. SUBSEQUENTLY GIVEN BY USER. 

//GO. OUTFIL DD DSN*A529E5 . VPI 17 . FI L El , 

// DISP*!NEW.PASS). 

// UMIT*TAPE. 

// LABEL r !19.SL. EXPDT *88 309), 

// VOL*! .RETAIN, . ,SER *200109), 

// DCB = {RECFM*VB,LRECL*38<4<4.BLKSIZE*16000, DEN *9) 

//GO.DD1 DD DSN*88G0SET,DISP*(0LD. DELETE) 

//■ THIS BLOCK MAY BE REPEATED FOR EACH FILE DESIRED WITH THE STEP • 
//■ CHANGED AND THE FILE ATTRIBUTES ADJUSTED. 

//STEP0002 EXEC PLIXCLO 
//PLI.SYSIN DD ■ 

« PROCESS SOURCE, FLAGU), GOSTMT j 
READTP . 

PROCEDURE OPTIONS! MAIN ) ; 

DECLARE SYSPRINT OUTPUT PRINT FILE ; 

DECLARE SYSIN INPUT RECORD FILE ; 


//GO. SYSIN 
// 

// 

// 

// 

// 
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Program TAPE TAPE continued 


TAPE 


TAPE 


DECLARE OUTFIL OUTPUT RECORD FILE » 

DECLARE DATA( 1920 ) FIXED BINARY! 15 ) l 
DECLARE COUNT FIXED BINARY! 31 ) > 

DECLARE £ND - 0F_FILE FIXED BINARY! 31 ) INITIAL! 0 
ON ENDFILE! SYSIN ) END_OF_FILE = 1 t 
ON ERROR SNAP BEGIN t 
ON ERROR SYSTEM t 
PUT DATA > 

STOP ; 

END ; 

DO WHILE! END_OF_FILE 
READ FILE! SYSIN ) 

WRITE FILE! OUTFIL 
END t 

END READTP t 


) I 


3 0 > » 

INTO! DATA ) < 

) FROM! DATA ) t 


/* 

//GO. SYSIN 
// 

// 

// 

// 

// 


DO DSN=FILE2# 

DISP»0LD, 

UNIT»TAPE, 

LABEL 3 !?, NL ) , 

VOL s SER=VPI17 , 

DCB=(RECFM=F, LRECL=38AO» BLKSIZE=38A0 > 


//GO. OUTFIL DD DSN=A529E5 . VPI17 . FILE2, 

// DISP 3 ! NEW, PASS), 

// UNIT=TAPE, 

// LABEL 3 ! 20, SL * EXPDT-S830A) , 

// VOL 3 ! , RETAIN* , ,REF=* . STEP0001 .GO. OUTFIL), 

// DCB 3 ! RECFM=VB.LRECL=38AA, BLKSIZE 3 16D00, DEN 3 A) 

//GO.DD1 DD DSN 3 88G0SET , DISP 3 ! OLD, DELETE) 

//STEP0003 EXEC PLIXCLG 
//PLI. SYSIN DD X 

* PROCESS SOURCE, FLAG! I), GOSTMT i 
READTP i 

PROCEDURE OPTIONS! MAIN ) > 

DECLARE SYSPRINT OUTPUT PRINT FILE i 
DECLARE SYSIN INPUT RECORD FILE i 
DECLARE OUTFIL OUTPUT RECORD FILE ; 

DECLARE DATA! 1920 ) FIXED BINARY! IS ) f 
DECLARE COUNT FIXED BINARY! 31 ) ; 

DECLARE END_OF_FILE FIXED BINARY! 31 ) INITIAL! 0 ) I 
ON ENDFILE! SYSIN ) END_.OF_FILE » 1 I 
ON ERROR SNAP BEGIN ; 

ON ERROR SYSTEM l 
PUT DATA t 
STOP > 

END > 

DO WHILE! END_OF_FILE 
READ FILE! SYSIN ) 

WRITE FILE! OUTFIL 
END i 

END READTP ; 

/* 

//GO. SYSIN DD DSN=FILE3, 

// DISP=OLD, 

// UNIT 3 TAPE, 


3 0 ) ; 

INTO! DATA ) 
) FROM! DATA 


) t 


28 


Program, TAPE TAPE continued 


TAPE 


TAPE 


LABEL»(3,Nl). 

V0L=SER=VPI17, 

DCB=! RECFM°F , LRECL a 38A0 . BLKSIZE s 38A0) 


// 

// 

// 

//GO. OUTFIL DO DSN=A529E5 . VPI17 .FILE3, 


DISP=!NEW,PASS>, 

UHIT=TAPE, 

LABEL 2 ! 21 , SL , EXPDT=8830A ) , 

VOL=! .RETAIN, , , SER*200109 ) , 

DCB=!RECFM=VB» LRECL a 38AA, BLKSIZE a 160QQ, DEN=A) 
DD DSN=88G0SET,DISP=! OLD, DELETE) 


// 

// 

// 

// 

// 

//GO.DD1 

//STEPOOQA EXEC PLIXCLG 
//PLI. SYSIN DD * 

X PROCESS SOURCE, FLAG(I), GOSTMT f 
READTP i 

PROCEDURE OPTIONS! MAIN ) ; 

DECLARE SYSPRINT OUTPUT PRINT FILE < 

DECLARE SYSIN INPUT RECORD FILE ; 

DECLARE OUTFIL OUTPUT RECORD FILE t 
DECLARE DATA( 1920 ) FIXED BINARY! 15 ) I 
DECLARE COUNT FIXED BINARY! 31 ) ; 

DECLARE END_OF_FILE FIXED BINARY! 31 ) INITIAL! 0 
ON ENDFILE! SYSIN ) EHD_QF_FILE ■ 1 i 
ON ERROR SNAP BEGIN ; 

ON ERROR SYSTEM > 

PUT DATA ; 

STOP > 

END t 

DO WHILE! ENDJJF^FILE 8 0 ) J 

READ FILE! SYSIN ) INTO! DATA ) t 
WRITE FILE! OUTFIL ) FROM! DATA ) j 
END ; 

END READTP ; 

/x 

DD DSN=FILEA, 

DISP=OLD, 

UNIT a TAPE, 

LABEL 8 ! A, NL ) , 

VOL e SER°VPI17 , 

DCB-!RECFM=F, LRECL-38A0, BLKSIZE = 38A0) 


) l 


//GO. SYSIN 
// 

// 

// 

// 

// 

//GO. OUTFIL DD DSN=A529E5.VPI17 . FILEA, 


// 

// 

// 

// 

// 

//GO.DD1 

// 


DISP=! NEW, PASS), 

UNIT=TAPE, 

LABEL 2 ! 22, SI , EXPDT e 8P30A) , 

VOL=!, RETAIN,,, REFsX.STEPOOOl. GO. OUTFIL), 
DCB : !RECFM-VB, LRECL°38AA, BLKSIZE=16000, DEN S A) 
DD DSN=88G0SET , DISP"! OLD, DELETE) 
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3. Program TAPE LISTE - read from the 3ystem tape and print the 
first 20 elements of every line. 


TAPE LISTE 


//B0589TL0 JOB 529E5, »TIM L SEN' . REGION 3 15 J6K, TIME 3 ! , 59) 

/JfJOBPARM LINES a 40 
/^PRIORITY IDLE 
//STEP0001 EXEC FORTVCO 
//FORT.SYSIN DD * 

INTEGER*2 DATA( 1920 ) 
i 10 CONTINUE 

READ! 5. END = 20 ) DATA 

WRITE! 6 , 99 ) < DATA! I ), I 3 1, 20 ) 

3 99 FORMAT! • 2016 ) 

GOTO 10 

J 20 CONTINUE 

STOP 

END 

1 /» DSN SHOWS NAME AND LABEL SHOWS POSITION OF FILE DESIRED, 

il //GO . FT05F001 DD DSN=A529E5.VPI12.FILEl, 

DISP=OLD, 

UNIT=TAPE» 

LABEL = (7 » SL ) < 

VOL=! » RETAIN# » »SER=200109) SERIAL • OF SYSTEM TAPE. 


// 

// 

// 

// 

// 


II t III I I I I I I I I I I I I I 

S? ^ i *5 2 ^ 2 5 2 2 N ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ — ^ 
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I - W * M «*• I If I»M i N I MN I m 

4«rm<fM*9lPO^KI9(W09aOtt*1«rviOMiilM«0>9t»9P<7N«f1 
_ ** ^ •'#*.••*»*.* a. *«n*r* 0 o-*M»<>***»*.*in«tfioM** 


I I I I I 1 I I I I I m I I I I • 1 I I 7 | I rnrntmm l m 1 | TT I l I 
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4. Program TAPE LISTL — read from the system tape and show the 
first few lines of the file in totality. 


TAPE 


LISTL 


//B0589TL0 JOB 529E5,'TIM OLSEN' , REGION 3 1536K,TIME S ( , 59) 

/XJOBPARM LINES=A0 

/^PRIORITY IDLE 

//STEP0001 EXEC FORTVCG 

//FORT.SYSIN DD * 

INTEGERK2 DATA< 1920 ) 

DO 10 I s 1 » 3 
READC5) DATA 
WRITE<6,99) DATA 

99 FORMATdH ,516/, 191(1016/), 516) 

10 CONTINUE 
STOP 
END 

/* 

//QO.FT05F001 DD DSN=A529E5.VPI12.FILE1, 

// 

// 

// 

// 

/✓ 


DISP»0LD, 

UNIT»TAPE, 

LABEL=(7 , SL ) , 

VOl=( , RETAIN, , ,SER=200109) 




Program TAPE LISTL continued - 


output 


Remaining 
line length 


File// 


i 


Day Minute 


Dummy 


5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 

5996 


Tm 

7560 ■ 
7660 * 
7700 * 
7750 • 
7810 • 
7870 
7930 
7990 
8050 
8100 
8150 
8210 
8260 
8360 
8390 
8650 
8500 


0 207 ’ id . 


Hour 


207 
•5856 -859 

‘6966 -671 

-6038 -277 

-8271 -1006 
-9016 -665 

-6985 -565 

-5809 -120 

-3235 -573 

-5266 -582 

-6599 186 

•5659 -139 


■6075 

■7003 


765 

793 


-6259 -231 

-5928 138 

-5533 1116 

-6779 599 


59 ' 
-1225 
-1309 
-6158 
-3573 
-2599 
-1756 
-1685 
-1392 
-19 
-1513 
-567 
-1060 
-3016 
-2293 
-866 
1011 
1610 


•5070 

•6268 

•6683 

■6972 

■5602 

■5015 

3990 

5255 

3131 

■3731 

-5685 

6581 

6366 

6529 

3223 

2272 

3617 


- 651 - 

- 1032 - 

- 2269 - 

- 1887 * 

> 1868 - 

- 1738 - 

- 670 - 

-19 

381 

232 

- 38 - 

- 558 - 

- 1622 - 

■ 1078 - 

- 930 - 

176 

817 - 


12975 

11726 

13269 

13159 

11781 

11909 

10255 

-9787 

-9539 

-9621 

11560 

11376 

10825 

12167 

10320 

-9180 

10173 


55 

-19 

-28 

-75 

0 

-56 

•112 

-19 

27 

-28 

-66 

18 

27 

66 

95 

27 

18 


606 

622 

622 

613 

613 

622 

622 

606 

622 

622 

622 

631 

631 

660 

622 

622 

622 


Time 

(in 

seconds) 


Time Channel 2 
1 


Once Square 
per wave 
Rev 


0 (5 zeroes indicate end of line) 


One line of original data (1920 elements) 
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5. Program POS CHECK - check number of data points in each cycle. 


POS CHECK 


//B0589TL0 JOB 5Z9E5.*T1- 'LSEN* ,REGI0N»1536K,TIME=(,59) 

/xLONGKEY HINGFAL 

/xJOBPARM LINES 3 AO. CARDS 1 100 

/xPRIORITY STANDARD 

//STEP0Q01 EXEC FORTVCG 

//FORT.SYSIN DO X 

CXXX 

CXXX POS CHECK READS SYSTEM TAPE* CHECKS POSITION VS SQ HAVE. 
Cxxx I, J. K COUNTERS; TEST 8 ON FOR CONTROL* POS 3 ROTOR 
C*xx ANGULAR POSITION. OPR = ONCE PER REV, SQH » SQUARE HAVE 
CXKK 

INTEGERX2 DATA( 1920 ) 

INTEGER I. J* K* POS, ON* TEST 
REAL SQH! 95500) 

CXKK 

Cxxx READ DATA FROM TAPE, ASSIGN SQUARE HAVE 

rxxx 

DO 10 1*1 * 500 
READ! 5 ) DATA 
DO 20 J s l, 191 

SQH((J-1)X191+J) » DATAUOXJ+5) 

20 CONTINUE 
10 CONTINUE 
Cxxx 

CXXX HEADING 
CXXX 

WRITE! 6 * 98 ) 

C 98 FORMAT! 1H1, AC POS SQHAVE’ ) ) 

98 FORMATCIHI,’ K POS SHCK) SH<K+1>»> 

Cxxx 

Cxxx CHECK SQUARE HAVE* ASSIGN ROTOR ANGULAR POSITIONS 

Cxxx HRITE POSITIONS AND CHANNEL VALUES 

CXXX 

POS 3 0 

ON = 0 

TEST = 0 

DO 30 K=1 * 95500 

IF (TEST. EQ. 1) THEN 
TEST = 0 
POS = 0 
ON 3 1 
ENDIF 

IF (SQH(K) .GE.6000) THEN 
IF ( SQH(K+1 ) . LE. 1000) THEN 
TEST 3 1 

HRITEC6.9Z) K, CPOS+1), SQH(K), SQH(K+1) 

92 FORMAT! 1H ,216, 2F7.0) 

ENDIF 

ENDIF 

IF (ON. EQ, 1 ) THEN 
POS 3 POS+1 

C HRITE!6,99) K, POS, TIME(K) , CH1!K), CHZ(K), CH3(K), 

C + CHA!K) , CH5(K), CH6(K), OPR(K), SQH(K) 

ENDIF 

C 99 FORMAT! 1H ,215, F8 . 0,8F7 .0) 

30 CONTINUE 
STOP 
END 


/X 

//GO. FT05F001 
// 

// 

// 

// 

// 


DD DSN=A529E5 . VPI10 . FILE1, 
DISP=OLD, 

UNIT=TAPE, 

LABEL = (l.SL), 

VOL=(, RETAIN, , ,SER=209109) 


34 


6, Program CYCLE SELECT - discard cycles where Che number of points 
is not equal to 360. 


CYCLE SELECT 


//BOS89TLO JOB 529E5, 'TIM OLSEN’ , REGION 8 1536K, TIME 3 ( * 59) 

/KLONGKEY HINGFAL 

/XJOBPARM LINES 8 AO >CARDS° 100 

/^PRIORITY STANDARD 

//STEP0001 EXEC FORTVCG 

//FORT . SYS1N DD * 

CXtfM 

CYCLE SELECT READS SYSTEM TAPE, FINDS 
I, J, K, L, M COUNTERS; TEST 8 ON FOR 
ANGULAR POSITION, OPR 8 ONCE PER REV, 

CYCLE 3 CYCLE NUMBER (IF ACCEPTED) 


CXXX 

Cxxx 

CXXX 

CXXX 

cxxx 


CYCLES WITH 360 UNITS. 
CONTROL, POS 8 ROTOR 
SQM 3 SQUARE WAVE, 


I NT EGER x 2 DATAC 1920 ) 
INTEGER I, J, K, L, M, 


Cxxx 

Cxxx 

cxxx 


REAL AP0S( 1146 ) , CH1C1146) 
REAL CH4( 1146 ) , CH5C1146), 


MNEH. POS, CYCLE, N 


CH2U146), CH3( 1146 ) 

CH6( 1146 ) , 0PRU146), SQWC1146) 


READ DATA FROM TAPE, ASSIGN CHANNELS 


N 8 720 
DO 10 1*1,6 
READ( 5 ) DATA 
DO 20 J°l, 191 
CHI ( ( I-l JX191+J ) 
CH2( ( 1-1 JX191+J ) 
CH3( ( I-l JX191+J ) 
CH4( ( I-l ) X191+J ) 
CH5<n-l)X191+J) 
CK6< (I-l )X191+J) 
0PR((I-1)X191+J) 
SQH((I-1)X191+J) 
CONTINUE 
CONTINUE 


DATAU0XJ-2) 
DATA(IOXJ-I) 
DATAC10XJ) 
DATA( 10XJ+1 ) 
DATAU0XJ42) 
DATA( 10XJ+3) 
DATA( 10XJ+4) 
DATA( 10XJ+5) 


20 
10 
Cxxx 

Cxxx HEADING 
Cxxx 


MRITE(6,98) 

98 FORMAT! ’1 CYCLE POS ORIG APOS CHAN1 CHAN2 CHANS ’, 
+ • CHAN4 CHAN5 CHAN6 1/REV SQWAVE’) 

CXXX 

CXXX CHECK SQUARE NAVE, ASSIGN ROTOR ANGULAR POSITIONS, 
Cxxx CHECK CYCLE FOR 360 POSITIONS, WRITE VALUES 
Cxxx 


POS = 0 

CYCLE 3 0 

DO 30 K=l»(N+360) 

POS 8 POS ♦ 1 
IF (SQW(K) .GE.6000) THEN 
IF (SQW(K+1).LE.1000) THEN 
IF (POS. EQ. 360) THEN 
CYCLE 8 CYCLE + 1 
POS 8 0 
DO 40 L=l,360 
M 8 K - 360 + L 
MNEW 8 (CYCLE-1 )X360 + L 
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Program CYCLE SELECT continued 


CYCLE SELECT 


99 

40 


+ 


APOS(MNEW) s MHEW 
CHl(MNEN) o CHI (M) 

CHZ(MNEW) a CH2CM) 

CH3CMHEH) a CH3CM) 

CH4<MNEW) » CH4CM) 

CHS(MHEN) s CH5<M) 

CH6(MNEN) a CH6CM) 

OPRCMNEW) a OPR(M) 

SQH(MNEW) a SQN(M) 

WRITE(6»99) CYCLE, L, M, APOSCMNEW), CHI (MNEW) , 

CH2(MNEW>, CH3CMNEW) , CH4CMNEW), CHS(MNEW) 
CH6( MNEW) , OPR(MNEW), SQH(MNEH) 
,I3,2IS,F6.0,7F7.0,F6.0) 


POS 
EHDIF 
ENDIF 
ENDIF 

30 CONTINUE 
STOP 
END 


FORMATOH 
CONTINUE 
ELSE 

0 


/# 

//GO. FT05FOOI DO DSN=A529E5. VPI10.FILE1, 
// DISP=OLD, 

// UNIT=TAPE, 

// LABEL = (1,SU, 

// VOL s < , RETAIN, , , SER s 200109> 

// 


# 
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. Program VOLT PLOT - plot data from each channel (1-6' 
wave, and once-per-rev signals against position 


square 


VOLT 


PLOT 


//B0589TLO JOB 329E5,'TIM OLSEN', RE0I0N»15S6K,TIME"(, 39) 

/XLONQKEY WIHQFAL 
/XJOBPARM LINES 3 40, CARDS 3 100 
/^PRIORITY STANDARD 
//STEP0001 EXEC FORTVCOV 
//FORT.SYSIN DD * 

CXMM 

C**x VOLT PLOT READS SYSTEM TAPE, PLOTS CYCLES WITH 360 UNITS. 

CM* I, J, K, L, COUNTERS) N 3 TOTAL • OF CYCLES, POS 3 ROTOR 
CXXN ANGULAR POSITION, OPR 3 ONCE PER REV, SQH = SQUARE HAVE, 

CXXX M 3 ORIGINAL DATA COUNT, MNEW 3 DATA COUNT OF ACCEPTED CYCLES, 
CX** CYCT » CYCLE COUNTER (IF ACCEPTED) 

CXMX 

INTEGERX2 DATA( 1920 ) 

INTEGER I, J, K, L, M, MNEH, POS, CYCT, N 
REAL APOS( 1146 ) , CHK1146), CH2( 1146) , CH3(1146) 

REAL CH4U146), CH5CU46), CH6(1146), OPRC1146), SQHC1146) 

N 3 3 

C**H 

C**» READ DATA FROM TAPE, ASSIGN CHANNELS 
CKX* 

DO 10 1=1, ( INK 360*(N+1 )/19l)+l ) 

READ( S ) DATA 


DATAC10KJ-2) 
DATA( 10XJ-1 ) 
DATA(IOXJ) 
DATAC 13*J+1 ) 
DATA( 10XJ+2) 
DATA( 1CXJ+S) 
DATA(10*J+4) 
DATA( 10KJ+5) 


DO 20 J»l,191 
CHI ( ( I-l )*191+J ) 

CH2( (I-l )*191+J) 

CH3((I-1)*191+J) 

CH4< ( I-l )*191+J ) 

CH5( ( I-l )X191+J ) 

CH6( ( I-l )*191+J ) 

0PR((I-1)*191+J) 

SQW( ( I-l )*191+J > 

20 CONTINUE 
10 CONTINUE 
CKKX 

C*KK HEADING 
C**K 

C HRITEC6 , 98) 

C 98 FORMAT< * 1CYCLE POS ORIG APOS CHANl 
C + ' CHAN4 CHANS CHAN6 1/REV 

CXXX 

C*XK CHECK SQUARE HAVE, ASSIGN ROTOR ANGULAR POSITIONS, 
CXXX CHECK CYCLE FOR 360 POSITIONS, HRITE VALUES 
CXXX 

POS 3 0 

CYCT 3 0 

DO 30 K 3 1, (<N+1 )X360) 

POS 3 POS + 1 
IF (SQH(K) . GE.6000) THEN 
IF (SQW(K+1) . LE.1000) THEN 
IF (POS. EQ. 360) THEN 
CYCT 3 CYCT + 1 
POS 3 0 
DO 40 L 3 1 , 360 
M 3 K - 360 + L 


CHAN2 CHANS 
SQHAVE* ) 
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Program VOLT FLOT continued 


VOLT PLOT 


MNEW a (CYCT-1 )X360 ♦ L 
APOS(MMEW) « MNEW 
CHI (MNEW) □ CHI (M) 

CH2CMNEN) » CH2(M) 

CH3CMMEH) = CH3(M) 

CH4(MNEW) » CH4(M) 

CHS (MNEW) a CHS(M) 

CH61MHEW) • CH6(M) 

OPR (MNEW) a OPR(M) 

SQW(MNEW) a SQW(M) 

C WRITE<6 ,99) CYCT, L, M, APOS(MNEW), CHI (MNEW) « 

C ♦ CH2(MNEW), CH3(MNEW), CH4(MNEW), CHS(MNEW), 

C ♦ CH6(MNEW) * OPR(MNEW) * SQW(MNEW) 

C 99 F0RMAK1H , 13,215, F6 . 0,7F7 . 0, F6 . 0) 

40 CONTINUE 

ELSE 
POS » o 
ENDIF 
ENOIF 
ENDIF 

30 CONTINUE 

CXXX 

PLOT SQUARE WAVE 


CXXX 

CXXX 


A XX « 

XFIR 

XINC 

AYY * 

YFIR 

YINC 


4.0 X N 
a 0.0 

a 90.0 

7.0 
a 0.0 

1000.0 


CXXK 

CXXK 

CXXX 


APOSCN+l ) a 0.0 
AP0S(N+2) a 90.0 
SQW(N+1) a o.O 
SQW(N+2) a 1000.0 
I TYPE = 0 
ITEXT = 11 
CALL PLQTS(0,0,50) 

CALL PLOT( 1 . , 1 . , **3) 

CALL VPISYM( 0.75,7.0,0.12, 'SQUARE WAVE*, 0. ,11) 

CALL AXIS( 0.0,0. OPPOSITION*, -8, AXX, 0.0, XFIR, XINC) 
CALL AXIS( 0.0, 0.0 .'VOLTAGE', 7, AYY, 90.0, YFIR, YINC) 
CALL LINE(APOS, SQW,N, 1,ITYPE, ITEXT) 

CALL PLOT(6 . ,6 . ,-999) 

PLOT ONCE PER REV 

AXX = 4.0 X N 
XFIR a 0.0 
XINC a 90.0 
AYY » 6.0 
YFIR = -15000.0 
YINC a 5000.0 
AP0S(N+1) = o.O 
AP0S(N+2) = 90.0 
OPR(N+l) a -15000.0 


38 



Program VOLT PLOT continued 


VOLT PLOT 


OPRCN+2) * 5000.0 
ITYPE » 0 
ITEXT » 11 
CALL PLOTS(0»0, SO) 

CALL PLOTC1 . >1 . , -3) 

CALL VPISYMC 0,75*6.0.0.12. 'ONCE PER REV',0.,12) 
CALL AXISCO. 0,0.0, 'POSITION* , "8 , AXX, 0 . 0, XFIR, XINC) 
CALL AXISC 0.0, 0.0. 'VOLTAGE', 7, AYY, 90. 0, YFIR, YINC) 
CALL LINE! APOS, OPR. N, 1, ITYPE, ITEXT) 

CALL PLOT(6 . ,6 . , -999) 

CXMM 

CHHX PLOT CHANNEL 1 
CXXK 

AXX ■•4.0 * N 
XPIR a 0.0 
XINC 9 90.0 
AYY 9 4.0 
YFIR 9 -16000.0 
YINC = 4000.0 
APOSCN+l ) = 0.0 
APOSC N+2) 9 90.0 
CH1CN+1 ) 9 -16000.0 
CH1CN+2) 9 4000.0 
I TYPE 9 0 
ITEXT = 11 
CALL PlCTSC0,0,5O> 

CALL PLOTCl.,1.,-3) 

CALL VPISYMCQ.75,6.0,0.'2, 'CHANNEL l',0.,9) 

CALL AXISC O.O,0.O,»POSr T ON', -8, AXX. 0.0, XFIR, XINC) 
CALL AXISC 0. 0. 0 . 0* 'VOLTAc 'j* ,7 , AYY, 90 . 0, YFIR, YINC) 
CALL LINE! APOS, CHI, N,l, ITYPE, ITEXT) 

CALL PLOTC6 . ,6 . ,-999) 

Cxxx 

C**M PLOT CHANNEL 2 
CXXN 

AXX 9 4.0 * N 
XFIR 9 0.0 
XINC 9 90.0 
AYY 9 5.0 
YFIR 9 -5000.0 
YINC 9 2000.0 
APOSCN+l ) 9 0.0 
APOSC N+2) 9 90.0 
CH2CN+1) 9 -5000.0 
CH2C N+2) 9 2000.0 
I TYPE 9 0 
ITEXT 9 11 
CALL PLOTS! 0, 0, 50) 

CALL PLOTCl.,1.,-3) 

CALL VPISYMC 0.75,6.0,0.12, 'CHANNEL 2*,0.,9) 

CALL AXISCO. 0,0.0, 'POSITION*, -8, AXX, 0.0, XFIR, XINC) 
CALL AXISCO. 0,0.0, 'VOLTAGE', 7, AYY, 90.0, YFIR, YINC) 
CALL LINECAP0S,CH2,N,1,ITYPE, ITEXT) 

CALL PL0TC6 . , 6 . , -999) 

CXXK 
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[j 

[1 

1 * 

• i Program VOLT PLOT continued 


VOLT PLOT 


Cxxx 

CXKX 


i 

1 


Cxxx 

cxxx 

cxxx 


1 

)! 

.M 


Cxxx 

cxxx 

Cxxx 


PLOT CHANNEL 3 

AXX a 4.0 X N 
XFIR a 0.0 
XINC = 90.0 
AYY a 6.0 
YFIR = -7000.0 
YINC a 2000.0 
APOSOI+l) = o.O 
APOS(N*Z) a 90.0 
CH3UI+1) a -7000.0 
CH3CN+2) a 2000.0 
I TYPE a 0 
ITEXT = ll 
CALL PLOTS<0,0,50) 

CALL PL0T( 1 . , 1 . » -3) 

CALL VPISYMt 0.75,6.0,0.12, 'CHANNEL 3',0.,9) 

CALL AXIS<0. 0,0.0. 'POSITION 1 , -8* AXX, 0 . 0 . XFIR, XINC) 
CALL AXIS<0. 0,0.0# 'VOLTAGE' ,7 , AYY, 90 . 0, YFIR, YINC) 
CALL LINE(AP0S,CH3,N, 1, ITYPE, ITEXT) 

CALL PL0T<6 . ,6 . , -999 ) 

PLOT CHANNEL 4 

AXX » 4.0 X N 
XFIR « 0.0 
XINC a 90.0 
AYY = 7.0 
YFIR « -14000.0 
YINC ° 2000.0 
APOS(NM) a 0.0 
APOSCN+2) a 90.0 
CH4<N+1) = -14000.0 
CH4( N+2) a 2000.0 
ITYPE « 0 
ITEXT » 11 
CALL PLGTSVQ, 0,50) 

CALL PLOT ( 1 . , 1 . , -3) 

CALL VPISYMCO .75,6 .0,0.12, 'CHANNEL 4',0.,9) 

CALL AXISC 0. 0,0.0, 'POSITION*, -fl, AXX, 0.0, XFIR, XINC) 
CALL AXIS<0. 0,0.0, 'VOLTAGE', 7, AYY, 90. 0, YFIR, YINC) 
CALL LINE(AP0S,CH4,N,1, ITYPE, ITEXT) 

CALL PLOT(6 . ,6 . ,-999) 

PLOT CHANNEL 5 

AXX a 4.0 X N 
XFIR « 0.0 
XINC a 90.0 
AYY = 5.0 
YFIR = -5000.0 
YINC = 2000.0 
APOS(Nn) a o.O 
AP0S(N+2) a 90.0 
CH51N+1) a -5000.0 


Program VOLT PLOT continued 


VOLT PLOT 


CH5IN+2) o 2000,0 
I TYPE « 0 
ITEXT « 11 
CALL PLOTSC0,0,50) 

CALL PLQT( 1. , 1. , -J) 

CALL VPISYM(0.75,6.0,0.12* 'CHANNEL 3', I 
CALL AXIS! 0 . Of 0 . Of •POSIT ION’ , -8.AXX, 0.1 
CALL AXIS( 0. 0, 0 . 0» * VOLTAGE* *7 , AYY,90 . 0 
CALL LINE! APOS, CHS, N, 1, ITYPE, ITEXT) 
CALL PL0T(6.,6.,-999) 

CXMH 

Cxxx PLOT CHANNEL 6 

CXXM 

AXX » A.O * N 
XFIR 8 o.o 
XINC s 90.0 
AYY = 6.0 
YF1R = -18000.0 
YINC 8 2000.0 
APOSCH+1 ) = 0.0 
AP0S(N+2) » 90.0 
CH6(N+1) a -18000.0 
CH6CN+2) a 2000,0 
ITYPE a o 
1TEXT a n 
CALL PLOTStO, 0.50) 

CALL PLOT! 1 . ,1 . »-3) 

CALL VPISYM 0.75, 6. 0,0. 12, 'CHANNEL 6', 
CALL AXIS(0. 0,0.0, 'POSITION' , -8, AXX, 0 . 
CALL AXISCO. 0,0.0, 'VOLTAGE' , 7, AYY, 90. 0 
CALL LINE< APOS,CH6,N,l, ITYPE, ITEXT) 
CALL PL0TC6 . ,6 . ,999) 

STOP 

END 

/X 

//GO . FT05F001 DD DSN=A529E5 . VPI10 . FILE1 , 

// DISPaQLD, 

// UNIT=TAPE, 

// LABEL =( 1 , SL ) , 

// V0L=(, RETAIN,,, SER=20Q109) 

// 
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. ,9) 

, XFIR, XINC) 
YFXR,YINC) 


l.,9) 

.XFIR, XINC) 
YFIR,YINC) 


